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Abstract 

$—1 1 

bX) This paper is devoted to study Bianchi type I cosmological model 

in Brans-Dicke theory with self- interacting potential by using perfect, 
anisotropic and magnetized anisotropic fluids. We assume that the ex- 
pansion scalar is proportional to the shear scalar and also take power 
law ansatz for scalar field. The physical behavior of the resulting mod- 
els are discussed through different parameters. We conclude that in 
contrary to the universe model, the anisotropic fluid approaches to 
isotropy at later times in all cases which is consistent with observa- 
tional data. 
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1 Introduction 



Many astronomical experiments and recent cosmological observations [T] in- 
dicate accelerated expansion of our universe. This expansion is believed by 
dark energy (DE), a cryptic exotic matter having large negative pressure 
that violates the strong energy condition. In radiation dominated era, the 
nucleosynthesis scenario elicits the decelerated expansion of the universe in its 
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early phase. To understand the nature of DE, many cosmological models like 
Chaplygin gas, phantom, quintessence and cosmological constant etc. have 
been proposed [2J. The modified theories of gravity like f{R) gravity, Gauss- 
Bonnet theory, higher dimensional theories of gravity, scalar tensor theories 
etc. have also been suggested [3J. Brans-Dicke (BD) theory of gravity is 
one of the most attractive scalar tensor theories due to its vast cosmological 
implications [1]. The varying gravitational constant acts as gravitational 
constant), the non-minimal coupling between the scalar field and geometry, 
compatibility with weak equivalence principle, Mach's principle and Dirac's 
large number hypothesis are some dominant features of this theory [5l [6]. 
The BD parameter should be constrained u > 40, 000 for its consistency 
with the solar system bounds [7]. 

Spatially homogeneous and anisotropic Bianchi type I (BI) model is used 
to study the possible effects of anisotropy in the early universe [8]. Some 
people [9] have constructed cosmological models by using anisotropic fluid 
and BI universe. Recently, this model has been studied in the presence of 
binary mixture of the perfect fluid and the DE [10]. Sharif and Kausar [TT] 
have discussed dynamics of the universe with anisotropic fluid and Bianchi 
models in f{R) gravity. Some exact BI solutions have also been investigated 
in this modifled theory |12] . 

In this paper, we construct solutions of the fleld equations for BI universe 
model in the presence of different fluids. The paper is organized as follows. 
In the next section, we formulate the fleld equations of BD theory for BI 
universe and some general parameters. Section 3 provides solution to the 
fleld equations in the presence of perfect fluid and then anisotropic fluid. 
The BI cosmological model with magnetized anisotropic fluid is investigated 
in section 4. A special case, m = 1, of the magnetized anisotropic fluid is 
also discussed. Finally, we summarize the results in the last section. 

2 Bianchi Type I Field Equations and Some 
General Parameters 

The BD theory with self-interacting potential is described by the action [13] 
^= /c/^Xv^[0i?-^0'°0,«-f/(0) + L^], a = 0,1, 2, 3, (1) 
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where loq and represent the constant BD parameter and the matter part 
of the Lagrangian respectively. Here we have taken SitGq = c = 1. Using the 
principle of least action, we obtain the field equations 

□0 = — ^ ^ ' . (3) 

3 + 2a;o 3 + 2u;o 

Here T^^,, T, □, A'^, represent energy-momentum tensor, its trace, box 
or d'Alembertian operator (□ = A'^A^), covariant derivative and the self- 
interacting potential respectively. Equation ([3]) represents the Klein Gordon 
equation or the wave equation for the scalar field. This theory reduces to 
general relativity (GR) when the scalar field is constant and the BD parame- 
ter is very large, i.e., w — oo However this is not true in general, e.g, the 
case of exact solutions. It is argued that this theory goes over to GR only for 
the non- vanishing trace of the energy- momentum tensor |l15]. For different 
values of this theory corresponds to other alternative theories of gravity. 
For example, it corresponds to Palatini metric /(-R) gravity, the metric /(i?) 
gravity and low energy string theory action for bj = —3/2, u = [16] and 
u = —1 [IT] respectively. 

The BI universe model is given by 



ds^ = - A\t)dx^ - B\t){dy^ + dz^), (4) 

where A and B are the scale factors. This model has one transverse direction 
X and two equivalent longitudinal directions y and z. The field equations ^ 
and ([3]) for the model (jlj) can be written as 

2AB fi2 ^ T^^^^_A B^ im 
AB^ B^ ~ 0+202 20 ' ^ ^ 

2^ + ^ = -Ill-^t^2^^-^ + ^ (B) 
B B^ 2 02 5 20 ' ^ ^ 

B A AB 2 02 ^A B^(P 20 

and the wave equation is 

^^r^^9^u ^ 2[/(0) - 0g 
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The corresponding average scale factor a{t), volume V and the mean Hubble 
parameter H are 



a{t) = {AB'Y/\ V = a\t) = AB\ H{t) =^-{^ + 2^). 
The directional Hubble parameters in x, y and z directions are given by 

Hx = ^, Hy = H^ = —. (9) 

The anisotropy parameter of expansion A and the deceleration parameter q 
are 

^^\ti^f, (10) 

1=1 

The isotropic expansion of the universe can be obtained for A = 0. The 
expansion and shear scalar turn out to be 

Since the field equations are highly non-linear, we assume power law for 
the scalar field (f){t) = 0o-B", a > for the expanding universe. For a spa- 
tially homogeneous metric, the normal congruence to homogeneous expansion 
implies that ^ is constant, i.e., "the expansion scalar B is proportional to 
shear scalar a" pjj. This leads to A = S™, m ^ 1 for BI model [18j, [20]. 
It is worthwhile to mention here that any universe model becomes isotropic 
when t — )■ +00, A — )■ 0, V ^ +oo, p > for the diagonal energy-momentum 
tensor. 



3 Anisotropic Fluid Model 

In this section, we first explore the BI model with the energy-momentum 
tensor of perfect fluid given by 

Tj^ = diag[p, -up, -up, -up], (12) 
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where p and uj represent the energy density and equation of state (EoS) pa- 
rameter respectively. Using this energy-momentum tensor in the field equa- 
tions ([2]) and ([3]), it follows that 

2?- + ^ = -^-^t^2^^-^ + ^ (U) 
B B"^ 2 02 B(f) (f) 0' 

where we have used A = B"^. 

The energy conservation equation for such a fiuid is 

p + 3H{l + oo)p = 0. (17) 
For < UJ < 1, this equation yields 

p = po5-('"+2)(i+-). (18) 

Subtracting Eq. f|T4|) from f|T5l) and using = 0o-B" (a > 0), we obtain 

5(t) = [(a + m + 2)(M + ^2)]'/^"+'"+'^; ^7^1 (19) 
where ki and A;2 are constants of integration. Consequently, we have 

A{t) = [{a + m + 2) {k^t + fc2)]™/('^+™+2) _ (20) 

Thus the model turns out to be 

ds'' = dt'- [(a + m + 2)(fcit + fc2)]'"'/^°+"'+'^rfa;2 

- [{a + m + 2){kit + k2)f^''^"'^^\dy^ + dz^). (21) 
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The corresponding parameters become 



u TT TT B mki 

^ B {a + m + 2){kit + k2) 



H 

e 



,m + 2 ki 



''{a + m + 2){kit + k2)- 
[m + 2)ki 



{a + m + 2){kit + k2) 
(m-l)2 kj 



3 '{a + m + 2y{kit + k2y 



(m + 2) 

Y ^ ^ [(^ ^ ^ ^ 2) {kit + fcs)] (^^+^ , 

_ d 1 3a + 2m + 4 

^ ~ Jt^W~ (m + 2) ■ 

Since a, m > (m 7^ 1), we have g > which yields the decelerated 
expansion of the universe. The mean anisotropic parameter of expansion 
(A = ^(^,^"2^^ ) is constant. 

In order to investigate the accelerated expansion model of the universe, 
we take the generalization of the perfect fluid, i.e., anisotropic fluid given by 

T"^ = diag[p, -p^, -py, -p^], (22) 

where p represents the energy density of the fluid while px, Py and pz denote 
pressures in x, ?/ and z directions respectively. Equation of state for this fluid 
is taken as, p = up, where EoS parameter uj may not be constant. By taking 
the directional EoS parameters uj^ = u) + 5, cu^ = a; + 7 and a;^ = w + 7 on 
x, y and z axes respectively, Eq.f l22|) can be written as 

T; = diag[l, -{u + 5), -(a; + 7), -(a; + 7)]p, (23) 

where 5 denotes deviation from a; on a; axis while 7 denotes deviations on 
y and z axis. Equation fl23|) with 5 = = 7 corresponds to the energy- 
momentum tensor for isotropic fluid. The energy conservation equation for 
the anisotropic fluid yields 

P + {l + uj){^ + 2|)p(t) + (5^ + 27|)p(t) = 0. (24) 
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By decomposing the anisotropic fluid into deviation free and anisotropy parts, 
we take anisotropy part equal to zero 



(5^ + 27|)pW = 0. (25) 

Since p 7^ 0, this implies that either both the deviation parameters S{t) and 
7 (if:) vanish or ^ = — For a more general solution, we take dimensionless 
deviation parameters as follows [21] 



, 2nB,A BA nA,A BA , , 

where n is a real dimensionless constant which describes the deviation from 
EoS parameter. 

The field equations for such fluid will be 

2?-^?^ = (^ + ^)P ^0'/'' ^(0) 

2 02 500 0' 



(28) 



Using Eqs.d^l]), and (121]) along with = 0o5", it follows that 
B , ,B^ n(m + 2)2^2 

^ + ^"^ + ^ + "^:B^- 3i^('^-^2)(^-i)0o = Q- 

Integrating twice, we obtain 

/CI, n(m + 2)2s~"^ 
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where and are integration constants. For B = T, x = X, y = Y and 
z = Z, Bl model turns out to be 

Some physical parameters are 

2fm-l)2 



(m + 2)2 



i/^ = mi7^ = m[2 .j.^im+i+a) 

^ ^ {m + 2) n/r-" 

3 Sa^ol?^ — 1) 

= 3H = im + 2)[2--— ^]T-(-+i+-), 

^2 ^ 2n/T-" 



3 30:00 (m — 1^ 

3 3 n/(m + 1 + 2a)T-('"+2+2«) 

{m + 2y ~ {m + 2y 3a<po{m - 1) 

- 2(m + 1 + a)T~(™+2+«))(2 ^^_iy(^+2+a)_ 

3q;0o("^ ~ 1) 

Since a, m > (m 7^ 1), these parameters except the deceleration pa- 
rameter, increase with the decrease in T and approach to zero as T — )■ 00. 
Also, for earlier times, the volume of the universe is zero while the expansion 
and shear scalar turn out to be infinite. For later times, the volume goes 
to infinite value while the expansion and shear scalar decrease to zero. This 
indicates that the universe expands from zero volume at infinite rate of ex- 
pansion. Since the anisotropy parameter of expansion is constant (it vanishes 
for m = 1), therefore the model does not isotropize for later times. In this 
case, the deceleration parameter q is found to be a dynamical quantity and 
can be negative for the appropriate values of the constant parameters. For 
later times and m > 1, the deceleration parameter turns out to be negative. 

The self-interacting potential U can be written from Eq. fl27p as follows 

^ U{T) = 20o((a + 2)m ^ + 1 + 2m)e^ ' s^oC— d-t- ) 

X T-(2™+4+") _ 2p. (32) 
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Equations and lead to 

dp 

co = -l ^ (33) 

{m + 2)p§ 

Substituting Eqs. flS^ and f l5^ in fl5U]) . we obtain 

p(T) = [0oa(m + 2)((3 + 2u;o)a(a + "^ + l)+4(l + 2m)+4a(m + 2) 

0oa'(m + 2)(3 + 2a;o) ^ 8a(m + 2) 

- 2a;o« ) — — — — + a-2 

z (3« — 2{m + zj) 

X [2T-("+2m+4)^g^^^ + 2) - (a + 2m + 4) (3a - 2(m + 2)))-^ 

- 2n/T-(2"+2™+4) (3a0o(^ - 1) (8a(m + 2) - (a + 2m + 4) (3a 

- 2(m + 2))))-] + 4c.(^ + 2)M^-l) ^^-(4+2.+2H(8^(^ + 2) 

3 

- (2a + 2m + 4) (3a - 2(m + 2)))"^ - n/T-(^+=^"+^")(3a0o(m - 1) 
X (8a(m + 2) - (3a + 2m + 4)(3a - 2(m + 2))))"^] + a^o 

^ (3 + 2a;o) (m + 2) (,^,^^4) n/r-(^"+^-+^) 
(3a-2(m + 2))^ 3a(/)o(m - 1) ^ 

+ ClT (3— 2(-+2)) , (34) 

where Ci is an integration constant. Inserting this value in Eq. fl52|) . one can 
obtain the corresponding self-interacting potential. The skewness parameters 
are given by 

2n(m + 2)(2 - Mp+^lflT )t-2{™+2+«) 

6{T) = — ^M^-^) ^ (35) 

3p 

-nm(m + 2)(2- 2n(m+2)^r-° xy-2(m+2+a) 
7(T) = ^ ^^oaim^l) > _ (3g) 

The deviation free EoS parameter (l33l) can be written as 

w(T) = -1- , . [[0oa(m + 2)((3 + 2a;o)a(a + m + l) 
p(m + 2) 

2n / o , 8a(m + 2) 



+ 4(l + 2m)+4a(m + 2)-2a;oa )-(a-2 



(3a-2(m + 2)) 

2(m + 2)(3 + 2u;o)i 



-][-2(a + 2m + 4)T-("+2m+5)(g^(^ ^ 2) 



2 

- (a + 2m + 4)(3a - 2(m + 2)))~^ + 2n/(2a + 2m + 4) 
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Figure 1: Plots represent energy density p versus time T for a < "^^"^^"^^ and 
a > ^'-"l"^^'' respectively. Here uq = 1.9, n = 2 and a = 1. 



X T-(2'"+2'"+^)(3a0o(^ - l)(8a(m + 2) - (a + 2m + 4) 

X (3a-2(m + 2))))-^] + L ^[-(4 + 2a + 2m) 

X T~(^+2"+2™)(8a(m + 2) - (2a + 2m + 4) (3a - 2(m + 2)))"^ + nl 

X (4 + 3a + 2m)r-(^+3"+2»n)(3^0^(^^ _ i)(8a(m + 2) - (3a + 2m 

+ 4)(3a - 2(m + 2))))-^] + "Vo|^^^^||^[-(« + 2m + 4) 



X 



y_(„+2™+5) ^ n/(2a + 2m + 4)T-(^"+^"-+5) 

3a0o("^ ~ 1) 



8a 1 8a 

^^(3a-2(m + 2))^ ^' ^ ^ 

where p is given by Eq. flM]) . The anisotropic expansion measure of anisotropic 
fluid is 

l_ _ '_ 3(/)oa(m-l) ' ^gg^ 

u 3ci;(T) 

Now we discuss the results for a < and a > H(!Z|l:Hl_ Figure 1 

indicates that the energy density is positive. For later times, it decreases 
and goes to zero for a > while it increases and approaches to infinity 

after big bang for a < ^iULt^ xhe self-interacting potential is positive only 
for a > as shown in Figure 2 and goes to zero for later times. Figures 

3 and 4 show that the skewness parameters 5(T) and 7(T) turn out to be 
finite at T = and approach to zero in future evolution of the universe for 
both cases. The anisotropy measure of expansion for anisotropic fluid goes 
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Figure 2: The self-interacting potential U{T) versus time T for « < ^^"^"^^^ 
and a > ^lll^t^l. Here ojq = —1.9, /3 = 2,n = —2, a = 1. 
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Figure 3: The skewness parameter 5{T) for a < '^^^^'^^ and a > 
respectively. Here uq = —1.9, a = 1, n = 2. 
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Figure 4: The skewness parameter 7(T) for a < '^^"^^'^^ and a > ^^"^g^^-* 
respectively. 
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Figure 5: The anisotropic measure of expansion parameter (5 — 7)/(X' for 
a < '^^''^^'^'> and a > '^'^"'■^'^^ respectively. 

to zero for T — )■ oo which shows that the anisotropic fluid approaches to 
isotropy for future evolution of the universe as shown in Figure 5. Notice 
that all these parameters decrease more rapidly with increasing values of the 
parameter m. 

At the initial epoch with a < ^^""^^^ and a > ^izili^)^ obtain u = —1 — 

and w = -1 + (^+2)(3a-2(m+2)) respectively. These indicate that 
for m > 1 and a > 0, the universe may be in phantom region or quintessence 
region. For later times with a < Hizi^tH) have a; = — 1 + — „w„^"„, — 

° 3 ' (m+2){3a—2(m+2)) 

and for a > follows that 
oj = -1-7 — ^-—r[[(poa{m + 2){{3 + 2ujo)a{a + m + l) +4:{l + 2m) 
+ Mm + 2) - 2a;oa ) - (« - 2 + ^^^-^^^^^) ^ 

X (3 + 2a;o)][-2(a + 2m + 4)] + + 2m 

+ 4)] [[0oa(m + 2)((3 + 2bjQ)a{a + m + 1) + 4(1 + 2m) + 4a(m + 2) 

2 0oa^(m + 2)(3 + 2a;o) ^ 8a(m + 2) , 
- 2 ((3«-2(m + 2))+"-2)][2 

X (8a(m + 2) - (a + 2m + 4)(3a - 2(m + 2)))"^] 
. (3 + 2a;o)(m + 2) ^_, 
(3a; — 2(m + 2)) 

This also shows that the universe will be in quintessence region or phantom 
region for later times depending on the value of the BD parameter. Thus the 
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model represents accelerated expansion of the universe. 



4 Magnetized Anisotropic Fluid Model 



In this section, we explore solution of the field equations for magnetized 
anisotropic fluid. We take anisotropic fluid with magnetic field along z axis 
and assume that there is no electric field. In this case, the scale factor Aii) 
is perpendicular to magnetic field while -B(t) is along the field lines. The 
magnetized anisotropic fluid is 



where pb represents energy density of the magnetic field. Using EoS for 
pressures in x, y and z directions as in the anisotropic fluid, Eq. (!39|) can be 
written as 

= diag\p + p^, -(w + 5) + p^, -(w + 7) - Pb, -(w + 7) - Pb], (40) 

where 5 and 7 are given by Eq. fl26|) . For 5 = = 7, Eq. fj40|) corresponds 
to the energy-momentum tensor for the magnetized isotropic fluid while it 
reduces to the anisotropic fluid for pB = 0. For 5 = = 7 Pb = 0, it 
represents the isotropic fluid. 

The field equations ([2]) and ([3]) for the model (j4]) and the energy-momentum 
tensor (l40l) become 



= diag[p + pb, -Px + Pb, -Py - Pb, -Pz - Pb] 



(39) 





2 02 50 

f/(0) 



{uj + 6)p-pB WO0^ r,B(}> 




(41) 



+ 



20 ' 



(42) 



(m + 1)— + m 




(u + 7)p + pB UJo 0^ 



52 



2 02 




(43) 



0+(m + 2)-0 



p(l - 3a;) - p{6 + 27) 2f/ (0) - 



(44) 



{2uo + 3) {2uo + 3) 
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where we have used the condition A = B^. 

The energy conservation equation for the magnetized anisotropic fluid 
yields pB = ^ along with Eq.f l2^ . Here /3 > is an integration constant. 
Subtraction of Eq. p2|) from fH3|) leads to 

2- H- 2(,„ + ! + .)(-)- ^^^^^ _ ^I^J -) ^ - ^^^^ Z^^^^^^ . (45) 

Taking i? = f{B), this turns out to be 

df* '^u . ^ IX* -4/9 

^ + -[(m + l + a)- 3^^^^_^^^J / = ^^(^_,^^.^3 , (46) 

where /* = and / = (m + 2)^ is a positive constant. This is the first- 
order linear non-homogeneous differential equation with variable coefficients 

2nlB~°' 

whose integrating factor is i?2(m+i+a)g3<^ga(„_i) _ After some manipulation, 
the solution becomes 

_4/35-("+2) 4nl/3B-^^'^+^^ 



(j)o{a + 2m)(m — 1) 3m(m — + 2m) 

30o("^ — 1)« 30o("^— l)a 



where C2 is an integration constant. This can also be written as 

_4/35-{"+2) 4n//35-2("+i) 



0o(a -|- 2m)(m — 1) 3m(m — 1)202^0; -|- 2m) 



X 



30o(m — l)a 30o("^ — !)« 



By taking B = T, x = X, y = Y and z = Z and using Eq. (H7|) . BI spacetime 
turns out to be 

as 



(j)Q{a + 2m){m — 1) 3m(m — l)20Q(a -|- 2m) 

^ 30o(m-l)a 3(/>o(m- l)a^^ 

- T^'^dX^ -T^{dY^ + dZ^). (48) 
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Now we discuss some physical features of this model. Since at T = 0, the 
scale factors will be zero, the model shows point type singularity 
The corresponding mean and directional Hubble parameters are 

= mHy = m\ 



H 



(j)o{a + 2m) (m — 1) 3m{m — iy(f)l{a + 2m) 
o(pQ[m — l)a o(po[m — l)a 

(49) 

(m + 2) . -4/3T-("+^) Anll3T-^'^"+^^ 



3 0o(a + 2m)(m — 1) 3m(m — l)20Q(a; + 2m) 
o(pQ{m — l)a o(po{'m — l)a 

(50) 

Since a, m > (m 7^ 1), these parameters increase with the decrease in T 
and approach to zero as T — )■ cxd. Also, these parameters take infinitely large 
values at T = 0. The remaining parameters are given by 

B = 3/i = (m + 2)[- 



4>o{a + 2m)(m — 1) 3m(m — l)20Q(a; + 2m) 

X (1 - ) + c2T-(™+^+")(i - j^^' )]v^ (51) 

30o("^ — l)a 30o('Ti — l)a 



m 



1)2 -4/3T-(-+4) 4nipT-^^"+^^ 



3 0o('^ + 2m)(m — 1) 3m(m — 1)20^(^q; + 2m) 
X (1- ' ) + C2T-2(-+2+")(l- )], (52) 

_ 3 , 3 , _4/3T~("+4) 4nZ/3 



X 



X 



(m + 2)^ 2(m + 2) >o(" + 2m)(m - 1) 3m(m - 1) 



2 



()q(q; + 2m) 30o("^ ~ 1)*^ 3(m — \)a 

_4/35-{"+2) 4nZ/35-2("+i) 



0o(a + 2m)(m — 1) 3771(171 — l)20Q(a + 2m) 



3(po{m — l)a 3(po{m — l)a 
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4(a + 2)/35-(°+3) 8(a + l)n//35-(2°+3) 



0o(« + 2m)(m — 1) 3m{m — l)202^a; + 2m) 

2>(f)Q{m — l)a 3m(m — 1)202(^0; + 2m) 30o(w^— l)a' 
2c2(m+l + a)5-(2-+3+2")(i ^n/E"" . 



+ C25-(-+^+")(|^^^^--)). 



In this case, the volume of the universe and anisotropic parameter of ex- 
pansion turn out to be the same as in anisotropic case. For initial time, 
the expansion and shear scalar become infinite while for later times, these 
decrease to zero. The deceleration parameter turns out to be a dynamical 
quantity. It can be negative for appropriate values of the constant parame- 
ters e.g., it becomes a negative for later times with m > 1. Notice that the 
expansion scalar, shear scalar and Hubble parameters are decreased by the 
component of magnetic field. 

We solve Eqs. pT]) and (jH]) simultaneously to obtain density p and the 
self-interacting potential f/(0) ~ U(T). The density is 

p{T) = (/)o(l + 2m + a(m + 2)-^)[- "^/^^ ' 



2 ''(a -t- 2m)(m - 1)00 



4n//3 T-(°+4) 2nlT- 



3m(m — 1)2 (a + 2m)0Q 30o(^ ~" !)« 
+ c.T-(^"""""Hl - )]-;^-^, (53) 



where U{T) is 



f/(T) = [((2a;o + 3)(a + m + l)a2_4a(l + 2m + a(m + 2)-^^))2(m 



+ 2) + a^m + 2)(3 + 2u;o)( ^,^''i"!,^ - (« - 2)) - 6a(l + 2m 

2[m + 2) — Sa 



+ aim + 2) - ^) .,^^f + ][4/3T-^[(c^ + 2m)(m - l)0o(8a 
2 2(m + 2) — 3a 



(m + 2) + 4(2(m + 2)-3a))]-^- ^^""^ 



3m(a + 2m) (m — 1) 



2^2 
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X (-T-(°+^)[(8a(m + 2) + (a + 4)(2(m + 2) - 3a))]-' + |n/0o ^(m 
- l)-^T-(2"+4)(8c^2(^^ + 2) + 2a(a + 2)(2(m + 2) - 3a))-' - 02(80 
X (m + 2) - (a + 2m + 4)(2(m + 2) - 3Q;))-^r-(^+2m+a) ^ 

3 

X (m - l)-^Q;r-(^°^+^'"+^^(8Q;(m + 2) + 2(q; + m + 2)(2(m + 2) 
^^^^,1 ^ Q;^(2a;o + 3) ^ -4(m + 2)^7"^ 4(m + 2) 



(2(m + 2) - So;) ^ (a + 2m) (m - l)0o 3m(m - l)^ 

/3n/ (.+4) _ 2n/r-^("+^)^o ^ („+2.,+4) 

(a + 2m)^ 3a{m-l) ' ^ 

2ni0^ pa+2n.+4)i , Mm + 2)^(1 - m)a ^ 4/3T-("+^)</>o ^ 
-1) 3 



3Q;(m — 1) 3 (Q; + 2m)(m — 1) 

X (8Q;(m + 2) + (Q; + 4)(2(m + 2)-3Q;))"^- 



4/3n/0o^ 



3m(Q; + 2m) (m — 1)^ 
X (-T-(2"+4)(8a(m + 2) + (2a + 4)(2(m + 2) - 3a))-^ 
2n/T"(^"+^) 

+ — (8Q;(m + 2) + (3Q; + 4)(2(m + 2)-3Q;))-^ 

o<po(m — ijct 

+ C2(-T-(^+2"»+2")(8a(m + 2) + (2a + 2m + 4)(2(m + 2) - Sa))"^ 

2^7T'-(3o!+2m+4) 

+ ^7T7 7^ — + 2) + (3q; + 2m + 4)(2(m + 2) - 3a))-^) 

3(po{m — l)a 

- 8a/3(m + 2)r-^(8a(m + 2) + 4(2(m + 2) - 3a))-^ + 6a(2(m + 2) 

- 3a)_i[^r-^ - 8^a(m + 2)T-^(8a(m + 2) + 4(2(m + 2) - 3a))-^] 

(l + 2m + a(m + 2) ^)[- ^ 



(2(m + 2)-3a)^ ^ ^ 2 ^^(a + 2m)(m - 1) 

4^"-^ .(2— (2a+3) _ 271/7'^^^"+^^ ^ _^ ^^^j.-(2a+2m+3) 



3m(m — l)2(a + 2m) 30o("^ — l)a' 

2n/T-(''^"+2™+3) 8a(m+2) 

- TT— +C3TW)-3a, 54 

3a(m — 1)00 

where C3 is an integration constant. Figure 6 indicates that the energy density 

is positive and decreases after big bang but it increases and approaches to 
infinity for later times with a < ^("^+^) . xhe energy density is positive but 
decreases to zero for any positive value of the parameter satisfying a > "^^"^^"^^ 
as shown in Figure 6. At the initial epoch, there is infinite energy density in 
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both cases as shown in Figure 6. Figures 7 indicate that the self-interacting 
potential remains positive in both cases (a < and a > ), The 

corresponding skewness parameters turn out to be 



2n{m + 2) . _4/3T-("+^) Anl(3T-'^^'^+'^'^ 



X 



1{T) 



3p (« + 2m)(m — 1)00 3m(m — l)^(a; + 2m)0o 

(1 - Jf^~], ) + c,T-(^+™+")(l - Jf^~], )], (55) 
o(po[m — l)a o(po[m — l)a 

-nm{m + 2) . _4/3T-("+^) Anll3T-^^''+^^ 



X 



3p (a + 2m)(m — 1)00 3m(m — 1)^(0; + 2m)0Q 

:i - ) + c,T-^(^+™+")(i - ^. )], (56) 

30o(m — l)a; 30o("^ ^ 1)" 



where p is given by Eq.f l53|) . Figures 8 and 9 indicate that the deviation 
parameters become finite at T = 0. For later times, these parameters con- 
verge to zero in both cases. From Eqs.f l24p and (125|) . the deviation free EoS 
parameter u can be written as 

u^iT) = -1 - -r^^. (57) 

The anisotropy measure of anisotropic fiuid, for the model f H8|) takes 
the form 

5-7 n(m + 2)2 -4/3r-('^+^) 4n//3T-2("+2) 



3u;(T) (a + 2m)(m — 1)00 3m(m — l)^(a + 2m)0Q 

X (1- ) + c.T-^^^+-+")(l- )]■ (58) 

30o(m — 1)q; 30o(m — l)a; 

Its behavior is shown in Figure 10. For initial epoch, this is finite while 
it goes to zero for the future evolution of the universe in both cases. This 
indicates that the anisotropic fiuid approaches to isotropy for later times. 
When T — > and a < ^lll^t^l obtain a; = — 1 — ('^+^'"+3") which shows 

3 ' (m+2) 

that the universe model will be in phantom region at initial epoch. For 
2(m+2) obtain u = —1 + -, — , — nvrr which shows that the 

3 ' (m+2)(3a— 2(r?i+2)) 

universe model will be in quintessence region at initial epoch. For later times 
with a < ^'•"V"^-' , it follows that u = —1 + -, — „w„^"„. — showing that the 

3 ' (m+2)(3«— 2(m+2)) ° 

universe may be in quintessence region. When a > Hizi^tHi^ the EoS parameter 
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Figure 6: Plots represent the energy density p{T) versus time T for a < 
and a > respectively. Here coq = —1.9, (3 = 2, n = —2. Green, 

red and blue lines show the graphs for m = 2, 3, 4 respectively. 




Figure 7: Plots show the self-interacting potential U{T) versus time T for 
a < and a > respectively. 
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Figure 8: The deviation parameter 5{T) versus time T for a < ^^"l"^^'' and 



a > 



2{m+2) 



respectively. Here loq = —1.9, (3 = 2 and n = —2. 
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Figure 10: Anisotropic measure of expansion versus time T is shown for 
a < and a > '^^"^^'^^ respectively. 

depends on the component of magnetic field and BD parameter indicating 
that the universe will be in phantom or quintessence region for appropriate 
values of the constant parameters. Thus, in each case, the model shows the 
accelerated expansion of the universe. 

Now we investigate a special case when m = 1. The scale factors become 
A{t) = B{t) = a(t) and the model turns out to be the FRW universe model 

ds^ = dt^ - a{tY{dx^ + dy"^ + dz^). 



Equation (H5l) yields 

a(t) = V2{^t + c,f\ (59) 
where C4 is an integration constant. The expansion scalar G turns out to be 

e = 3iJ = 3 -) = 3W— (W— t + C4) . 
a V on V on 
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This shows that the Hubble parameter and the expansion scalar are constant 
at earlier time. As time increases, both of these parameters decrease indi- 
cating expanding universe in its earlier time. From Eqs. fHT]) and (jS]), the 
energy density p(t) and the self-interacting potential U(t) are 



Pit) 



U{t) 



So 



'2(3 ,c-4 



-[2(f)oa{{3 + 2uo)a{a + 2) - 12(1 + a) + 2uoa'^) + 



a\3 + 2uo) 
{2 -a) 



^'o(« + 2) 



160o«^ 
(2-«) 



(3(1 + a) 



woo^ 2£ 
2 '^^y 3n 



. (g. + 2)-' + + ^"°>^ - t°y^°':(3(i + a) 



3n{2 - a) 



3n{2 - a) 



Here C5 is an integration constant and the constants Si and S2 are given by 



^1 = -(l/2)[20oa((3 + 2a;o)a(a + 2) - 12(l + a) + 2a;o«2 

,,2/Q , n. , \ A IRA ^.2 , , ^2 



+ 



X 



a (3 + 2wo)0O/ , r,\2 1600^"' , ^ woa' 



(2-«) 



a + 2r- 



(2-«) 



(3(1 + a) - 



,2/3 



^2 



^fp(«+2)-^ 

3n 

2("-^)/Voa^(3 + 2a;o)/3 _ 4a0o/32"/^ 

^ 6n(2 - a) 6n(2 - a) 

u;oa2 2(-"2)/2^ 



(3(1 + «) 



+ 0o(3(l + a) 



3n 



Some other parameters are 



-2/3 /2^ 2 
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5 — 7 



+ 



+ 



+ 




Uj{t) 



where a ^ 2. We discuss two cases: < a < 2 and a > 2. Clearly, 
the energy density is constant at initial epoch and approaches to infinity 
for later times in both cases. The anisotropy parameters are constant at 
initial epoch and go to zero for later times. Likewise the anisotropy measure 
of expansion of anisotropic fluid ^ approaches to isotropy for later times. 
The anisotropy parameter of expansion is zero as m = 1. At the initial 
epoch, the deviation free EoS parameter shows that the universe may be in 
quintessence region by choosing appropriate values of the constants in both 
cases. For later times with < a < 2, we obtain uj{t) = —1 — '^^ich 
indicates that the universe will be in phantom region. For o; > 2, it follows 
that uj{t) — —1 — '■"g^-* , which also shows that the universe will be in phantom 
region for future evolution. 



In this paper, we have constructed the Bl universe models in BD theory of 
gravity with perfect, anisotropic and magnetized anisotropic fluids. We have 
constructed exact solutions in each case. For anisotropic and anisotropic 
magnetized fluid models, the physical behavior of the energy density, self- 
interacting potential, skewness parameters and anisotropy parameter of ex- 
pansion of anisotropy fluid have been plotted for non-zero value of n with 
a < "^^"^^"^^ and a > ^i!^|tHl. The results are summarized as follows. 

• In the case of anisotropic as well as magnetized anisotropic fluids, the 
skewness parameters and anisotropic measure of expansion of anisotropic 



5 Summary and Discussion 
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fluid go to zero indicating the isotropic behavior of the fluid for the fu- 
ture evolution of the universe. This result coincides with those already 
available in literature for Bianchi type III model in f{R) theory [TT] 
and Bianchi type iVI)o model in GR [23] . 

In each case, the energy density remains positive. All the figures in- 
dicate that energy density increases after big bang and approaches to 
infinity for later times with a < in both anisotropic as well as 

magnetized anisotropic fluids. For a > ^iindll^ it decreases and goes to 
zero in both cases. 

For anisotropic fluid, the self interacting potential is positive only for 
a > ^^"y"^-* and decreases to zero for later times while for magnetized 
anisotropic fluid, it remains positive in both cases. 

All the physical parameters H, H^., Hy, G and a increase with the de- 
crease in T and go to zero as T — )■ oo. These parameters take inflnitely 
large values at T = 0. In contrast to the perfect fluid, the decelera- 
tion parameter for anisotropic fluids is a dynamical quantity and can 
be negative for the appropriate choice of constant parameters, in par- 
ticular for later times with m > 1. This corresponds to accelerated 
expansion of the universe. 

In the anisotropic magnetized fluid, all the physical parameters are 
reduced by the component of magnetic fleld with n > 0. 

The deviation free EoS parameters indicate that the universe may be 
in quintessence or phantom region at initial epoch as well as for later 
times for an appropriate values of the constant parameters in all cases. 
Thus the models represent the accelerated expansion of the universe. 

The anisotropy parameter of expansion is constant (vanishes for m = 1) 
indicating the model does not isotropize for later times in all cases. 

A special case m = 1 for the magnetized anisotropic fluid has also been 
discussed which yields FRW universe model. In this case, the deviation 
free EoS parameter indicates that at initial epoch, the universe may 
be in quintessence region while for the future evolution, it will be in 
phantom region. 
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It would be interesting to construct exact solutions in the presence of anisotropic 
fluid for other Bianchi models in BD theory. 
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